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Kaon Weak Matrix Elements with Wilson Fermions* 
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We present results of several numerical studies with Wilson fermions relevant for kaon physics. We compute 
the Bk parameter by using two different methods and extrapolate to the continuum limit. Our preliminary result 
is _B^? S (2 GeV) = 0.66(7). AJ = 3/2 A' -t tytt matrix elements (MEs) are obtained by using the next-to-leading 
order (NLO) expressions derived in chiral perturbation theory (ChPT) in which the low energy constants (LECs) 
are determined by the lattice results computed at unphysical kinematics. From the simulation at /3=6.0 our (pre- 
liminary) results read: (7T7r|0^(2 GeV)]K)j =2 =0.14(l)(l) GeV 3 and (ivTY\of S (2 GeV)|K) /=2 =0.69(6)(6) GeV 3 . 



1. K°-K° MIXING 

The main problem of the computation of the 
K°-K° mixing amplitudes with Wilson fermions 
is the lattice operator mixing (not present in the 
continuum) introduced by the explicit breaking 
of the chiral symmetry due to the Wilson term. 
In order to compute the amplitudes one has to 
subtract contributions of operators with differ- 
ent naive chiralities. Two alternative propos- 
als which avoid the operator subtractions have 
been recently suggested: the first uses twisted 
mass QCD [1] while the second uses suitable 
Ward Identities (Wis) [2]. A first comparison of 
the standard method (with subtractions) to the 
Wis method (without subtractions) has been pre- 
sented in ref. [3] for the case of the parameter Bk, 
defined as 

(K a \ Ol \K )^lf KmK BK, 

with Oi =Qi-fii = s q 7m (1 - l5 )d a s b lti (l - 75 )d b . 
Qi is the parity conserving part of 0\ and un- 
der renormalization mixes with operators of dif- 
ferent naive chiralities, while Q\ is its parity vi- 
olating part and due to CVS symmetry is mul- 
tiplicatively renormalizable. The basic idea for 
the Wis method is based on the following WI ob- 
tained from a T3 axial rotation: 



m^(n°( 2 )P(G,q)Q 1 (O)P(t„q))+0( a ) = 

(P{t x ,(i)Q 1 {0)P(t y ,p)) + (rot. sources), 

where P(t, p) = Y^ x e~ lpx d(t, x)7 5 s(£, x) and 
Id = dj^d — U75U. SU(3) symmetry together 
with charge conjugation C imply that the terms 
involving the rotation of the sources are identi- 
cally zero. Instead of computing the three point 
function involving Q\ one can compute the four 
point function involving Qi, thus avoiding lattice 
operator mixing. From now on the procedure to 
obtain Bk is common to both methods and can be 
found in ref. [4]. Note, however, that 0(a) effects 
present in both methods may be different. Our 
simulation has been performed at three values of 
/3, by keeping roughly the same physical volume. 



p 


6.0 


6.2 


6.4 


volume 


1&* x 52 


24 a x 64 


32' 1 x 70 


# confs. 


500 


200 


135 


a" 1 from K/K* 


2.05(7) 


2.68(11) 


3.44(10) 


aM P 


0.30-0.39 


0.20-0.35 


0.15-0.23 



*Talk presented by Mauro Papinutto. 



The values of B^ 8 (2 GeV) for each value of P are 
given in Table 1, where the last column contains 
the values obtained from a linear extrapolation 
to the continuum limit (see Fig. 1). The average 
of the two gives Bf 5 ^ GeV) = 0.66(7). We are 
currently preparing a simulation at a larger (3 to 
discard, eventually, the point at (3 = 6.0 which 
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p 


6.0 


6.2 


6.4 


oo 


a (fm) 


0.096(3) 


0.073(3) 


0.057(2) 





no subt. 


0.818(32) 


0.763(25) 


0.750(37) 


0.634(95) 


stand. 


0.890(45) 


0.826(34) 


0.819(40) 


0.70(12) 



Table 1. Values ofB^ s (2 GeV) in function of a 
could lie outside the region linear in a. 




without subtraction 
• with subtraction 



weighted average: B K " 5Ba '(2 GeV] - 0.66(7) 



a (GeV 



Figure 1. Continuum limit for B^ s (2 GeV) 



2. K 



MEs IN THE CHIRAL LIMIT 



The K^tt MEs of 3/2 ,0 3 s /2 (defined in ref. [6]) 
can be related, in the chiral limit, to those of 
the left-right operators of the AS=2 basis as ex- 
plained in ref. [6] . By using the soft pion theorem 
(nn\O 7 , 8 \K ) I=2 = -i/(V2U)(Tr+\0 3 7 / J\K+), we 
obtain the A7=3/2 MEs of the electro- weak pen- 
guin (EWP) operators 7 ,0% (defined in Eq. (1) 
of ref. [5]) in the chiral limit. Our results (in the 
Ms scheme at 2 GeV) read: 



6.0 6.4 

{■kit\O 7 \K ) I=2 [GcV) 0.132(10) 0.136(14) 0.136(10) 
(7r7r|O 8 |jr )j = 2[GeV 3 | 0.658(44) 0.604(44) 0.533(34) 

For Os discretization errors seem quite large and 
we are currently investigating this feature. 

3. AI = 3/2 K -> tttt AMPLITUDES 

The procedure to extract these amplitudes 
by using ChPT at NLO has been explained 
in ref. [5,7,8]. We are interested in K — » tttt 
matrix elements of the EWPs and of the op- 
erator O4 (defined in Eq. (1) of ref. [5]) which 
mainly determines the physical K + — > 7r + 7r° 
decay amplitude. In practice we work with an 
unphysical kinematics (called SPQR) where the 
kaon is at rest, one pion is at rest and the other 
has either momentum or momentum 2tt/L. 
We have 3 pion masses and 3 kaon masses for 
each pion one (M w € [0.51,0.88] GeV, Mk € 



t 



• M fl (M K ,M t ,EJ p=0 
- M S (M H ,M„,E„) p=2n/L 



Figure 2. Numerical data for the ME of 8 

[0.51, 1.46] GeV), and the energy is in general 
not conserved (it can be injected through the 
weak operator [9]). We fit the numerical data 
to the expressions obtained at NLO in quenched 
ChPT [5,8]. For example, for the EWPs we have 



M s 7 



SPQR 
8 



Ml 



-IS.. 



X 7,S 



■(chiral logs) 



2 " 2 I ^■'• U --- L ^ 1 



(2M| + AMI) 



{AMI + ™l) 



(1) 



where 7 7,8 is the LEC of the LO while the 5j' 8 are 
the LECs of the NLO. For the moment we have 
not considered the finite volume correction (the 
Lellouch-Luscher (LL) factor [10]) which is not 
yet known for the SPQR kinematics. In this way 
we determine some combinations of the LECs of 
the NLO and we use them to compute the values 
of the physical amplitudes (see refs. [5,8] and refs. 
therein). Again in the case of the EWPs we have 



1 + 



{A-Kff 

f 7.8 



(chiral logs) 



sl' 8 ) 



+ 25^ + 28^ 

1 ^l7,8 . r7,8 . ,<c7,8 1 ,c7,8 . <-,r7,8\ 2 

+ ( S[ +o 2 + 4o 4 ' +45 5 ' + 25 6 ' J m^. 

For the EWPs, 6 counterterms of the NLO con- 
tribute in the SPQR kinematics and, as shown in 
Eq. (1), we can determine only 4 combinations 
of the corresponding LECs. These are enough to 
compute the physical amplitude. For 4 there are 
again 6 counterterms at the NLO but all of the 
corresponding LECs have to be separately deter- 
mined. Consequently, in the second case the fit is 
much more complicate. 

For the EWPs, the analysis at LO (which is 
0{p ) and thus correspond to fit the numerical 
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OS — full logarithms 
# 

i r 



0.5 ■ • M S (M K ,M.,EJ/FIT p=0 

©M^M^M^EJ/FIT p=2ir/L 



OS - quenched logarithms 



-1* *tt n *♦ 



0.5 • M B (M K ,M,,EJ/FIT p=0 

o M^M^M^EJ/FIT p=2it/L 



Figure 3. Fits to Ms with full and quenched logs 
data to a constant) gives the results 



M*r, .EL < 



1.5 GeV 1.0 GeV 0.8 GeV 



(itit\0 7 \K») i=2 

(7TTT\Os\K ) I = 2 



0.076(6) 
0.909(45) 



0.080(6) 
0.860(45) 



0.094(6) 
0.803(46) 



where we varied the number of data points used 
in the fit by choosing only points with masses and 
energies below a certain cut-off. The dependence 
of the results on the cut-off is due to the fact that 
higher orders in ChPT arc not negligible (as one 
can argue from Fig. 2 in the case of Og)- We thus 
introduce the NLO 0(p 2 ) contribution in the fit. 
In Fig. 3 we show the quality of the fits obtained 
by using cither the logarithms of full ChPT or 
those of quenched ChPT. We see that, according 
to what expected [11], only the quenched loga- 
rithms have the appropriate form to fit the data. 
On the other hand, the fit is not very sensitive 
to the quenching parameters S and a [12] and, 
in the lacking of their precise determination, we 
set their values to be 5 = 0.1 and a = 0.0. The 
contributions of the LO and NLO to the fit are 
shown in Fig. 4. The fit now shows much less de- 
pendence on the cut-off and we quote the results 
obtained with Mjf,£ T < 0.8 GeV as our central 
values, and use those with other cut-offs to esti- 
mate the systematic error due to the higher-order 
corrections in the chiral expansion. This leads to 
the preliminary results (in which the finite- volume 
LL corrections have not been included): 

(HOp(2 GeV)\K)i= 2 = 0.14(1)(1) GeV 3 , 

(7T7r|0;p(2 GeV)\K) I=2 = 0.69(G)(6) GeV 3 . (2) 

Note, however, that we can not claim to be sen- 
sitive to the chiral logarithms. In fact we could 
have fitted our numerical data by including only 
the counterterms of the NLO and not the log- 
arithms. Also in this case we would have ob- 
tained a very good fit (of quality comparable to 
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Figure 4. Various contributions to the fit of A4s 

that which includes the chiral logarithms) and re- 
sults which depend very weakly on the cut-off and 
which are very close to those obtained by includ- 
ing the logarithms: (0 7 )i=2 = 0.13(1) (1) GeV 3 
and (O 8 )/ =2 =0.64(6)(2) GeV 3 . We also note that 
this determination and the one in Eq. (2) are both 
in reasonable agreement with the results obtained 
from K-^tt MEs in Sec.2 (at /3=6.0). Finally, it 
is important to remark that the analysis at LO 
would lead to a large systematic error (as large as 
40% in the case of 7 ). 

Also in the case of O4 the NLO contribution 
turns out to be large. The fit is more problematic 
due to the large number of parameters and the 
errors on some of the NLO LECs (and therefore 
on the whole NLO contribution) is quite large. 
As an indicative result (obtained with Mk, E n < 
1 GeV) we quote .A(A+^7r+7r°)=0.0097(18) LO + 

0. 0039 .85) NLO GeV 3 (the experimental value is 
_4 cxp =0.0104 GeV 3 ). We are presently trying to 
consider other kinematical configurations which 
could improve the quality of the fit. 
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